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^ \ Abstract 

, ^ , '■ We employ perturbation analysis technique to study multi-asset 

\ portfolio optimisation with transaction cost. We allow for correlations 

CO ' in risky assets and obtain optimal trading methods for general utility 

\ functions. Our analytical results are supported by numerical simula- 

■ tions in the context of the Long Term Growth Model. 

•• 1 Introduction 

in ■ 

o ■ 

I In recent years, the study of portfolio optimisation under non-zero transac- 

■ tion cost has received its due attention (Davis &: Norman 1990; Atkinson & 
> ! Wilmott 1995; Morton & Pliska 1995; Akian et al. 1996; Atkins & Dyl 1997; 

; Atkinson & Al-Ah 1997; Atkinson et al. 1997; Atkinson k Mokkhavesa 2001, 

■ 2003, 2004; Mokkhavesa & Atkinson 2002; Chellathurai & Draviam 2005). 



In the hterature, it is found that the incorporation of transaction fees into 
the model introduces a no-transaction region around the original optimal 
curve, surrounded by purchase and sale regions (cf. figure [1]). Most of pre- 
vious work focused on having only one risky asset (stock) and one risk-free 
asset (bond), except the study by Atkinson and Mokkhavesa (2004) in which 
portfolio with multiple risky assets is analysed. In this work, the authors 
are able to obtain the optimal investment strategy with the assumption 



'Present address: Citigroup Global Markets Asia Limited, Citibank Plaza, 3 Garden 
Road, Central, Hong Kong 

^Author for correspondence (C.Leel@physics.ox.ac.uk). 



1 



Figure 1: A schematic diagram depicting the optimal trading strategy for the 
long term growth model. The optimal holding of the risky asset, ^*(cx 11), 
is shown by the red curve. The sale-no-transaction boundary is given by 
{A* + a+) and the purchase- no-transaction boundary is given by {A* -|- a_) 
where a_ < (cf. equation II ip . In this work, we assume that a+ = — a_. 




that the risky assets are uncorrelated. Here, we go beyond this restriction 
and consider correlated risky assets. By assuming that the purchase and 
sale boundaries are of an equal distance away from the optimal curve, we 
obtain analytical expressions for the optimal trading strategy for general 
utility functions. We further support our analytical results with numerical 
simulations in the context of the Long Term Growth Model. 

The plan of the paper is as follows: In §2, we will consider portfolio 
optimisation without transaction costs, thereby introduce the dynamic pro- 
gramming method employed. We will then consider trading with transaction 
cost in §3. The details of our simulation method in support of our analytical 
results are given in HI For reference, the expressions for the derivatives of 
the value function in terms of the expansion parameter are given in [H 

2 Trading without transaction costs 

We consider a market with investment opportunities on n stocks and a risk 
free bond, and we let Ai(t), B{t) and U{t) = B{t) + YA^^ ^{t) be the values 
held in stock i, the value held in risk free bond and the total wealth at time 
t respectively. We assume that Ai{t) follows a geometric Brownian motion 
with growth rate and volatility ai, and the risk free bonds, B, compounds 
continuously with risk free rate r. The volatilities <Tj, growth rates //j and 
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interest rate r are assumed to be constant. Cash generated or needed from 
the purchase or sale of stocks is immediately invested or withdrawn from the 
risk free bonds. In the absence of transaction costs, the problem is easily 
solved without recourse to perturbation analysis and this section will serve 
to familiarise the readers with the use of dynamic programming method in 
this optimisation problem. 

The market model equations are represented by the foUowings: 

dAi = fiiAidt + UiAidXi , z = 1, • • • , n 



dB 



= fiiAidt + UiAidXi , i 
= rBdt = r lu-^Ai 



dt 



dU = rBdt + J2 fJ'i^i^i + J2 f^iAdXi , (1) 

i=l 1=1 

where Xi , i = 1, . . . ,n, are Weiner processes whose correlations, — 1 < pij < 
1, are assumed constant. At time i = 0, an investor has an amount Il{t = 0) 
of resources and the problem is to allocate investments over the time horizon 
t G [0,r], so as to maximise the following expectation value: 



F{U{T)) + / I{U{t'))dt' 



E 



The functions / and F can represent anything from utility to the year end 
bonus of the trader. For example, if we assume that / = and F(Il(T)) = 
log(n(r)), then the opimisation problem constitutes the Long Term Growth 
Model and the goal would then be to optimise the logarithm of the final 
wealth. To make financial sense, we will assume that the utility functions 
are increasing and concave down, i.e., 

dl 



du 



> 



OF 



dip 

d^F 

dip 



< 



< . 



(2) 
(3) 



We restate the optimisation problem in dynamic programming form by 
first defining the optimal expected value function, J(n, t): 

F{U) + I^^I{U{t'))dt'\ . (4) 



J{Il,t) = maxEj 



We now apply the Bellman Principle and Ito's Lemma to the above value 
function to obtain the following Hamilton-Bellman-Jacobi equation (Kamien 
1991): 







max 




+ 2^t^iAi Qjj + 2 ^HjAiAj 
i=l i,.i=l 



(5) 
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with the boundary condition J(n, T) = F{Il(T)). In the above equation, 
Qij = GiGjPij is the standard covariance matrix, and Ai act as the control 
parameters in the context of dynamic programming. 
In matrix notation, we can rewrite equation ([5]) as 



= max 

A 



Jt + I + rUJn + {ft ■ A)Jn + ^(^ • ^^^)^nn 



(6) 



where symbols without an index denote the corresponding vectors (e.g. A = 
{Ai, . . . , An)^)- We have also introduced a new vector fi, which is defined 
to be (/ii - r, ...,//„ - r)^. 

By differentiating equation ([6]) with respect to A, one obtains as the 
solution to the HBJ Equation: 

dJ ^ d'^J^. ^ 

^/. + ^OA = 0. (7) 
Therefore, the optimal portfolio corresponds to: 

2.1 Example: the Long Term Growth Model 

In this model, our aim is to maximize E[logn(T)]. The value function is 
thus 

J(n,i) =maxEt[log(n(r))] . (9) 

such that J(n, T) = log(n(r)). This boundary condition together with the 
differential equation obtained by substituting equation ([H]) into equation ^ 
implies that 

j(n, t) = log n + (^r + ^/i • n-^p}j (T-t) . (lo) 

The optimal portfolio from equation ^ is therefore given by: 

A*=Un-^il. (11) 

For the case of having two-risky assets, the optimal portfolio corresponds to 
having B* = qU with q = I — Tr[r2~^/i], and A* = piU with pi = 
By the Ito's lemma, we have 

d(log n) = (^rq + fiipi + fi2P2 - + (3dX (12) 

where 



P = y O'lH + 2o-iCr2PiP2Pl2 + (T2P2 ■ (13) 



4 



Figure 2: Analytical and simulation results for the two-risky-asset market 
with model parameters given infH Left plot: Yi = (log(n(T))) denotes the 
performance of the optimal trading strategy, while Y2 denotes the perfor- 
mance of the sub-optimal strategy obtained with the correlation between 
the two risky assets ignored. Right plot: The plots of (Yi — Y2) and the stan- 
dard errors of the means for Yi and Y2 versus time based on our simulations. 
Note that S = 4000 is the number of samples in the simulations. 




The optimal expected payoff in this model is therefore: 

E[log(n(r))] = E[log(n(0))] + (^rq + /iipi + fi2P2 - ^ • (14) 

We now consider a two-risky-asset market model. With the model pa- 
rameters given inUl the optimal stock holdings in this case are Al = 0.067x11 
and A2 = 0.467 x 11 (cf. equation (jlip ). The performances of this optimal 
trading strategy based on our analytical expression in equation p4p and our 
numerical simulations (cf. H] for details of simulation method) are given in 
figure [21 If the correlation in the risky assets is ignored, the optimal portfolio 
becomes: = 0.3 x 11 and A2 = 0.5 x H. The corresponding performance 
is shown to be sub-optimal in figure [2j 



3 Trading with transaction cost 

We will now include transaction cost into our discussion. As the transaction 
cost usually amounts to a small percentage (~ 0.5%) of the total transaction, 
we employ perturbation method to analyse this optimisation problem with 
the transaction cost as the expansion parameter. By keeping track of the 
first few lowest order terms, we will derive the first order correction to the 
optimal trading strategy determined under no transaction cost. 
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We assume that the transaction fee is proportional to the asset under 
transaction and the proportionahty constant is denoted by k. Note that we 
again define the total wealth, 11, as 



U = B + Y.^i 



(15) 



i=l 



The market model equations in this case are: 

dB = rBdt-{l + k)dLi{t) + {l-k)dM,{t) 



dn 



U-{l-k)J2A dt - (1 + k)dU{t) + (1 - k)dMi{t) 

i=l 

HiAidt + dLi{t) - dMi{t) + aiAidXi , i = 1, • ■ • , n 

n 

rn + (1 - A;) E ( - rAidt + fI^A,dt + OiAidXi 



2 = 1 



11 



i=l 



(16) 



where Li{t) and Mi{t) represent the cumulative purchase and cumulative 
sale of assets Ai during the time interval [0,T]. The optimal expected value 
function J(n, A, t) is as before: 



J(n, A, t) = max E 

Li, Mi 



FiU{T))+l^^mit'))dt 



(17) 



and the corresponding HBJ equation is (Kamien 1991): 




I aj V- / . dLi dMA dJ 



+ 



i=l ^ 
n \ n 



dt dt J OA 



i=l ) i=l 



dt 



dt 



dJ_ 
dU 



" ^ , , /i av id^j d^j 



2dAidAj 2 5n2 dAidU 



(18) 



Here, L, and Mi are the control parameters from the dynamics programming 
perspective. 

3.1 Three regions 

By isolating terms involving dL or dM separately in equation (jlSp , we arrive 
at three separate cases: 



Casel: ^-kll<0^nd-^-kll>0. 
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In this case, the maximum in equation (llSp is achieved by choosing dLj = 
and dMj = oo, which is equivalent to selhng at maximum rate. 



Case ^;|^-fc|^>Oand-|^-fc|^<0. 

In this case, the maximum is achieved by choosing dLj = oo and dMj = 0, 
which is equivalent to buying at maximum rate. 

Case 3: 1^ - < and - fcf^ < 0. 

In this case, the maximum is achieved by choosing dLj = and dMj = 0, 
which indicates that no transactions are needed. 

We note that it is not possible to have ~ ^Wi ^^'^ ~'M' ~ ^§Ti 
be both greater than zero as we assume that J is an increasing function 
of n. This can be broadly interpreted as more wealth cannot decrease the 
value function from the trader's point of view. 

With the above consideration, the optimal trading strategy can be seen 
to be partitioned into three separate regions: sale, purchase and no-transaction 
regions (cf. figure[T]). In other words, if the portfolio is in the sale (purchase) 
region, the optimal strategy is to sell (buy) stocks until the portfolio is at 
the no-transaction region boundary, and thus bring the portfolio back into 
the no-transaction region. Inside the no-transaction region, dL and dM are 
identically zero and hence J satisfies the HBJ equation with k = 0. 

3.2 Continuity and optimality assumptions 

To make progress with our analysis, we will assume that the optimal value 
function, J, is everywhere continuous and that its derivatives are also con- 
tinuous. We call the latter the optimality assumption. The validities of 
these assumptions are discussed in Morton & Pliska (1996) and Whalley & 
Wilmott (1997). 

We now restrict ourselves to one risky asset for notational convenience. 
Suppose that the point (11, A, t) is inside the sale region, when a very small 
quantity of assets, h is sold, the risk-free bond increases by the amount 
h{l — k), while the whole portfolio value is reduced by kh. As /i — > 0, the 
value function J must be the same after the sale (the continuity assumption), 
we therefore have 

lim J {U + kh, A, t) = lim J{U, A - h,t) (19) 

^.^^^ y{U + kh,A,t)-J{n,A,t) ^ ^.^ J{U,A-h,t)-J{U,A,t) 
h^o kh h^o h 

k^ - (21) 
By a similar argument, we can conclude that inside the purchase region. 
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we have 



k- 



dJ dJ 



(22) 



dli dA ■ 

By applying again the same argument to equations (f2T]) and (f22]) with 
the use of the optimahty assumption, we have that in the sale region and at 
the sale-no-transaction boundary: 



5V 



-k- 



dJ 



(23) 



dAidn ' 

and in the purchase region and at the purchase-no-transaction boundary: 



dA?~' 



dAidU 



(24) 



Inside the no-transaction region, the value function, J, must satisfy equa- 
tion ([18]) with dL = dM = 0, i.e., 



i=l 



i=l 



i=l 



dJ_ 



" ^ , /l 9V lQ2j Q2j 



2dAidAj 2 5n2 Sylvan 



(25) 



These equalities are to be supplemented by the boundary condition at t = T: 

j(n,AT) = F(n). 



3.3 Perturbative expansion and order matching 

We now redefine the Ai coordinate as Ai = A*(n,t) + k^^^ai, where A* 
is the optimal value of stock i held when k tends to zero, and introduce 
the modified value function, H, such that H{Il,a,t) = J(n.,A,t). InHJ we 
display the various derivatives of J in terms of H and a. 
We further expand H(Il, a, t) in powers of k^^'^ as: 

Ho{IL,a,t) + k^'^ Hi{n,a,t) + k^'^ H2{Ti,a,t) 
+ kH3{U,a,t) + k^/^H4in,a,t) + 0{k^/^) . (26) 

The reason for expanding H and Ai in powers of k^^^ is out of necessity and 
has previously been studied in the literature (Atkinson & Wilmott 1995, 
Rogers 2004). 

We will from now on keep track of the expression up to the first non- 
trivial correction: 0{k^^^). By matching the orders of k, equations ()21|) and 
(|22p at the sale-no-transaction boundary (corresponds to the + sign in it) 
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and at the purchase-no-transaction boundary (corresponds to the — sign in 
±) become: 

= ^ , 0<m<2 (27) 
and equations ([23l) and (f24]) become: 



= , 0<m<2 (30) 



daf daidn y Sa^ J • ^""^^ 

Inside the no-transaction region, after expanding H according to equation 
()26p and cohecting terms of the same order in A;, we arrive at the following 
conditions: 

1. C(/c~2/3) Equation: VHq = 0, where V is an operator defined as 

Ei:,=iA,a2^,^. with 

1 8A* dA* " 1 r)A* " 

^^^■ = 2 an 9n ^ ^^M^^^r+2^^.A*^,*-^E^^'^A*^;; (33) 

h,l=l h=l 

2. 0(A;-i/3) Equation: Pi^i = 0. 

3. C(l) Equation: X'i72 = —MHq, where A4 is an operator defined as 

(n \ n 1 " 

i=l / i=l id=^ 

4. C(A:V3) Equation: VH3 = - ELi ai9A,(X-ffo) - MHi. 

5. 0(A:2/3) Equation: VH^ = w Sm=i " -^^2- 
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Combining the 0{k~'^^^) equation with equations (I27p and ()30p when 
m = 0, one finds that Hq is independent of a. Combining the 0{k~^/^) 
with equations ()27p and (I30p when m = 1 shows that Hi is independent of 
a. Combining the 0(1) equation with equations (j27p and (j30p when m = 2 
shows that H2 is independent of a. The 0{k^^^) equation together with 
equations (I28|) and (I3ip imply that H3 is independent of a. In summary, 
by matching the coefficients of the various orders in k, we determine that 
Hq, Hi, H2, H3 are independent of a. 

Without loss of generality, we focus on the first asset and let ai+ denotes 
the width of the purchase-no-transaction boundary, and ai_ the width of 
the sale- no-transaction boundary (cf. figure [T|) . From equations (j29p and 
(f30l) . we find that at the boundary A* + ai_: 

dHi dHo , , 

+ ^ = (35) 



9a 1 9n 

d'^H. 



dal 

and at boundary A* + ai+, we have: 

dH^ dHo 
dai dU 

d^H, 



daf 



(36) 

(37) 
(38) 



As we have established that Hq and H2 are independent of a, with the 
0{k'^^^) equation, we can conclude that H4 has the following general form: 



i?4(n,a,t) = ^hj{U,ai,t)ai , (39) 

j=0 

where aj denotes the set {om : m > 1}. In other words, H^ is a polynomial 
in ai with a degree of at most four. We now make the simplifying assumption 
that a+ = —a-. This is equivalent to saying that the transaction (buy or 
sell) boundaries are of the same distance away from the unperturbed optimal 
curve. We note that this assumption is proved to be true in the case of 
having uncorrelated risky assets (Atkinson and Mokkhavesa 2004). With 
the assumption of equal magnitude, we can conclude that /13 = at the 
boundaries by subtracting equation (f36l) from equation ([381) . In particular, 
we have 

6/i4a+ + /i2 = . (40) 

By summing equations ()35p and ()37p . we can further determine that hi = 
at the boundaries. By subtracting equation (j35]) from equation ([37|) . we 
conclude that ai+ satisfies: 

dHo 



dU 



4/1401 , + 2/i2ai+ . (41) 
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Substituting equation (fiOl) into equation (|^T]) . we obtain 

1 dH^ 



8/14 dli 



(42) 



To calculate h^, we invoke the 0{k'^/^) equation: By comparing the coeffi- 
cient of the a\ term on both sides, we find that: 



/14 



^ 2iDn ■ 
So finally, ai± can be expressed as: 



(43) 



a 



i± 



3Du dHo fd^Hp^ 
' al dll\ dip 



(44) 



where Hq is the optimal value function when transaction cost is absent. 

In general, denoting the trading boundary for stock i by aj+, we have 
the following general expression for the widths of the trading boundaries: 



a,; 



-1 



a? dli \ dip 



1/3 



(45) 



For any financial model where Hq is known, the above equation together 
with equation ^ provides an analytical description of the optimal trading 
strategy. This is the main result of this paper. 

3.4 Example: the Long Term Growth Model 

According to equations [10] and [HI 



logn+ (r + ]^fi.VL-^(^ (T-t) 



Ho{U,t) 

A* = iiii ^/i 
Combining these with equation ([33]) . we have 



(46) 
(47) 



1 dA* dA* 



n,iAlAt + -nijA*A* 



dA. 



* n 



2 dU dU ^ ^" ' ' 2 * ^ 5n , , 

h,l=l h=l 



J2 ^^hA*AU48) 



= tf|i(A-f)-V + ^n[f^~V]?-/ii[f^-V]f} • (49) 
The width of the boundary for stock i is therefore (cf. equation (fiSl) ): 

1/3 



n 



3k 



1 



-1 ,-.l2 



(50) 
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Figure 3: Simulation results for the two-risky-asset market with transaction 
cost. Left plot: Yi = (log(n(r))) denotes the performance of the optimal 
trading strategy, while Y2 denotes the performance of the sub-optimal strat- 
egy obtained if the correlation between the two risky assets is ignored in 
the calculations for the boundary widths. Right plot: The plots of (Yi — Y2) 
and the standard errors of the means versus time based on our simulations. 
Note that S = 15000 is the number of samples in the simulations. 
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If the risky assets are uncorrelated, the expression above coincides with the 
result of Atkinson & Mokkkhavesa (2004). 

We now employ this optimal trading strategy to the two-risky-asset mar- 
ket considered before. The optimal curve corresponds to: = 0.067 x 11 
and A2 = 0.467 x 11 (cf. equation (fTTIl ). and according to equation (f50l) . the 
boundaries widths are: «!+ = 0.167 x k^/^U and 02+ = 0.710 x k^/^U. The 
performance of this strategy is shown in figure [3l If we ignore the correlation 
between the risky assets in calculating the boundary widths, ai+ and 02+ 
become 0.508 x k^^^U and 0.760 x k^^^U respectively. The trading strategy 
employing these boundaries together with the same optimal curve as before 
is shown in figure [3] and can be seen to be sub-optimal, albeit the difference 
is small. 

In figure HI the portfolio's temporal evolution of a particular simulation 
is shown together with the transaction amounts displayed. 

4 Conclusion 

In conclusion, we have employed perturbation method to study multi-asset 
optimisation for arbitrary utility functions. By making the assumption that 
the sale and purchase boundaries are of the same distance away from the 
optimal curve, we arrived at an analytical expression for the optimal trad- 
ing strategy. We have also supported our analytical results with numerical 
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Figure 4: A particular simulation run with the optimal trading strategy in 
the two-risky asset model. Upper plot: the temporal evolution of values held 
in bond and stocks. Lower plot: The transactions performed according to 
the optimal trading strategy in the time interval t S [4,6]. The purchases 
(sales) of stock Ai are denoted by blue (green) crosses, and the purchases 
(sales) of stock A2 are denoted by black (red) triangles. Note that the 
transaction amount is not infinitesimal only because of the discrete time 
evolution (tick time) in the simulations (cf. U]). 
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simulations in the context of the long term growth model. 
Details of numerical simulation method 

We consider a portfolio consisting of two risky assets and one risk-free 
asset. The values in the bond and risky assets are updated as follows: 

B{t + At) = B{t) +rB{t)At (51) 
Ai{t + At) = Ai{t)+ niAi{t)At + aiAi{t)y/Ai zi{t) (52) 
A2{t + At) = A2{t) + H2A2{t)At + a2A2{t)^/Ai Z2{t) (53) 

where zi{t) and Z2{t) are random numbers drawn from the normal distribu- 
tion with zero mean and a standard deviation of one, such that the correla- 
tion coefficient between zi{t) and Z2{t) is pi2- 

In the case of trading without transaction costs, the portfolio is updated 
after each iteration according to equation (llip . When transaction costs are 
present, trading only occurs when the value of the risky assets are outside 
of the no-transaction region (cf. figured]), i.e., if 

A,{t) i [4*(t) +a,+ (t),4*(t) +a,_(t)] . (54) 

When such an event occur, the portfolio is adjusted such that Ai{t) is moved 
back to the nearest boundary and the cost of transaction is subtracted from 
the wealth. For example, if A\(t) > Al{t) + a+(t), then the portfolio is 
adjusted so that: 

B{t) ^ B{t) + {l-k)[Ai{t)-Al{t)-aMt)] (55) 
Ai{t) ^ A*(i) + ai+(t) . (56) 

The simulations always start with a total wealth of 1 at the optimal 
portfolio distribution and the set of parameters employed are: r = 1,/xi = 
1.3, H2 = 1-5, ai = (72 = 1, pi2 = 0.5, k = 0.005 and Ai = 5 x 10"^ 

Change of variables Letting H(U, a, t) = J(n, A, t) with A = A*{Il, t) + 
fc^/^a, we have the following expressions for the derivative of J in terms of 
H and a. 

1=1 ' 

(60) 



d^J _ d'^H " / d'^H dA* dH d'^A* 

dip ~ dYP~ daidU dU ^ 9a, dU^ 
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n 



d'^H dA*dA* 
daidaj dH dU 



(61) 



n 



daidaj dU 



(62) 
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